Assessing the Impact of Continuous Vaccination and Voluntary Isolation on the Dynamics of COVID-19: A Mathematical Optimal Control of SEIR Epidemic Model

In order to study the impact of continuous vaccination and voluntary isolation for the COVID-19, a susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model is proposed. A basic regeneration number R0 is defined to determine the extinction or persistence of the disease. We numerically analyze the impact of key parameters based on actual parameters of COVID-19, such as the vaccination rate, population importation rate, and natural (or causal) mortality transmission rate on the dynamics of disease transmission. Then we obtain sensitivity indices of some parameters on R0 by sensitivity analysis. Finally, the stability of the system and the effectiveness of the optimal control strategy are verified by numerical simulation.


Introduction
Since the outbreak of the coronavirus disease 2019 (COVID-19), it has been significantly impacting the world economy and our lives. Globally, as of 4:00 pm CET, 28 February 2022, there have been 434,154,739 confirmed cases of COVID-19, including 5,944,342 deaths, reported to the WHO. As of 26 February 2022, a total of 10,585,766,316 vaccine doses have been administered.
With the emergence of multiple mutant strains of novel coronavirus, including Alpha, Beta, Gamma, Delta, and Omicron, many scholars further explore the transmission and epidemic pattern of COVID-19 from different view. ese research results are helpful to learn about the transmission and mode of infectious diseases and provide reliable information for the prediction and control of infectious diseases.
In the study of the dynamics of infectious diseases, the susceptibility-infection-recovery (SIR) model was mainly established by Kermack and McKendrick in 1927 with the kinetic method. In order to better study the characteristics of infectious diseases and obtain the best prevention strategies, scholars propose many improved models based on the traditional SIR model for different situations. By introducing time variation factors, Chen et al. [1] present a new SIR model, which has stronger adaptability and robustness in predicting the number of confirmed cases and inflection points of infectious diseases. Wang et al. [2] propose a SIR model with time-varying isolation protocol. In this model, the dynamic evolution of infectious diseases can be comprehensively analyzed and evaluated. In addition, Hota et al. [3] introduce a closed-loop framework combined with the SIR model and expound the significance of early detection through two feasible optimization questions.
In view of the latent characteristics of some infectious diseases, SEIR models are further proposed based on the SIR model. Li et al. [4] fit the basic regeneration number curve based on the SEIR model. Trends in COVID-19 outbreaks in China, the USA, India, and Iran are predicted and analyzed.
e results show that the rapid and efficient isolation measures adopted in China are significant in suppressing COVID-19. In addition, some new models take into account the influence of different factors on infectious diseases, such as the presence of temperature and humidity [5,6]. It is concluded that COVID-19 is more infectious and lethal at low temperatures and humidity. Besides, many research studies show that the different epidemic prevention measures, such as mandatory isolation, wearing masks, vaccination, and government control policies, have significant effects on preventing the rapid spread of COVID-19 [7][8][9].
Due to the problem of asymptomatic infected persons in infectious diseases, Yu et al. [10] propose an SEIR-AQ model by considering factors including prevention and control efforts, isolation strategies, and asymptomatic infected persons. e SEIR-AQ model can well anticipate the spread trend of COVID-19, providing technical support for the scientific assessment of the infectious disease situation. Analysis from the SEIR-AQ model shows that prevention and control isolation, medical follow-up isolation, and other measures have a significant inhibitory impact on the spread of COVID-19. In addition, some scholars use deep learning methods to study the impact of the COVID-19 outbreak on human society [11]. Some other scholars have analyzed and studied the transmission trends of the COVID-19 outbreak. e spatiotemporal evolution and transmission trends of the COVID-19 epidemic are analyzed using an ontological modeling approach by Liu et al. [12].
Although control measures, vaccination, latent infectivity, and other factors are considered in these studies, the stability and existence of the model are not discussed from the mathematical view point. erefore, we establish a new susceptible-exposed-infected-recovered-quarantine-vaccines (SEIR-QV) model based on the latent period and the influence of prevention and control measures of COVID-19. In mathematical theory, the extinction or persistence of a disease can be determined based on the basic reproduction number and the Lyapunov function. We numerically simulate and analyze the impact of key real-world parameters of COVID-19. What's more, our sensitivity analysis indicates that the related parameters have significant effects on the stability and existence of the SEIR-QV model.

The SEIR-QV Model
In this paper, the susceptible-exposed-infectious-recoveredasymptomatic-quarantine-vaccines (SEIR-AQV) model is developed based on the traditional SEIR model. e SEIR-AQV model takes into account the effect of isolation based on the traditional model and divides the population into unisolated susceptible people (S), unisolated exposed people (E), unisolated infected people (I), isolated susceptible people (S q ), isolated exposed people (E q ), asymptomatic infected people (A), hospitalized people (H), and recovered people (R).
e SEIR-AQV model has the following assumptions: first assumption: medical resources are adequate, i.e., isolated exposed people can be directly converted to hospitalized people after diagnosis; second assumption: once infected with COVID-19, one must undergo inpatient treatment to recover, i.e., infected people are unlikely to recover on their own; third assumption: both isolated exposed people and hospitalized people are isolated from the outside world, i.e., they are not infectious; fourth assumption: the infected people will be immune after recovery, i.e., he or she will not become susceptible again after recovery. e bin transformation relationship is shown in Figure 1. e equations of the SEIR-AQV model are as follows: where c means the rate of exposure; η means the natural mortality rate; q represents the proportion of isolated; δ represents the vaccine coverage rate; β represents the probability of transmission; θ stands for the ratio of A relative to the transmission capacity of I; υ stands for the ratio of E relative to the transmission capacity of I; b 1 , b 2 , b 3 stand for the conversion of I, A, andE q , respectively to H; λ 1 represent the rate of conversion of S into S q ; λ 2 stands for the rate of conversion of S q into S; σ stands for the ratio of E to I conversion and σ � 1/Incubation period; e means the ratio of E to I conversion; α 1 , α 2 , α 3 mean the rate of conversion of I, A, and the rate of cause-specific death of H; r means the rate of recovery of H; and Λ means the rate of population importation. However, given the reality of the COVID-19 epidemic, the status of isolated susceptible people is unlikely to persist over time, and isolated susceptible people have no translational relationship to hamlets other than unisolated susceptible people. erefore, the population type of isolated susceptible people and the parameters involved in it are removed. Considering also that the proportion of asymptomatic infected people is very small and has minimal effect on the overall transmission trend of the COVID-19 outbreak, the population type of asymptomatic infected people and its involved parameters are also removed. erefore, we improved the SEIR-AQV model in the context of the COVID-19 outbreak. We construct a new seasonal susceptible-exposed-infected-removed-quarantine-vaccines (SEIR-QV) model with the population transformation relationship bin view shown in Figure 2.

2
Computational Intelligence and Neuroscience e equations of the SEIR-QV model are as follows: where α 1 , α 2 represent the rate of cause-specific death of unisolated infected I and hospitalized H, respectively, and b 1 , b 2 represent the rate of conversion of unisolated infected I and isolated exposed E q to hospitalized H, respectively.

Proof of Stability of the SEIR-QV Model.
To study the stability of equation (2), it is sufficient to study the stability of the following three formulas in (2) [13,14]: Considering the biological significance of the system, the dynamical properties of equation (3) are discussed mainly in the closed set Ω � (S, E, I)∈ R 3 + S| ≥ 0, E ≥ 0, I ≥ 0 , where R 3 + denotes the first trigonometric limit of R 3 and contains the boundary. When

Basic Regeneration Number.
e basic regeneration number (R 0 ) indicates the number of people infected by a patient during the average disease period when all are susceptible at the beginning of the disease. R 0 � 1 can be used as a threshold to decide whether the disease is extinguished or not. When R 0 < 1, the disease will become extinct. When R 0 > 1, the disease will persist. e basic regeneration number is closely related to the stability of the endemic equilibrium point.
Next, we study the basic regeneration number R 0 of equation (3). Let x � (E, I, S) T , then equation (3) can be rewritten as where e Jacobi matrices of F(x) and V(x) at the disease-free equilibrium P 0 � ((Λ/η), 0, 0) are, respectively, where e basic regeneration number, denoted by R 0 , is thus given by the following equation: 4 Computational Intelligence and Neuroscience
e other characteristic roots satisfy the following equation: It is known that So we can obtain It is easy to verify for R 0 < 1 and R 1 < 1 that the roots of the quadratic equation are as follows: All have negative real parts, i.e., the disease-free equilibrium point P 0 � ((Λ/η), 0, 0) of equation (3) is proved to be locally asymptotically stable.

Theorem 2. For equation (3), the disease-free equilibrium
Proof. By considering the Lyapunov function, Clearly, the solution of V(t) along equation ( Take en, where the inequality sign is obtained based on S ≤ S 0 . When R 0 ≤ 1, there is dV(t)/dt ≤ 0. us, dV(t)/dt ≤ 0 when and only when erefore, if R 0 ≤ 1, the disease-free equilibrium point P 0 of equation (3) is globally asymptotically stable in Ω. □ Theorem 3. S * , E * , I * are positive when R 0 > 1.

Proof
It is known that the parameters are positive and q < 1, so S * > 0 obviously holds the following equation: It is known that R 0 > 1, q < 1 and the parameters are positive, so E * > 0. 6 Computational Intelligence and Neuroscience It is known that R 0 > 1, q < 1 and the parameters are positive, so I * > 0.

Numerical Simulation.
e stability of disease-free and endemic equilibrium points is discussed in this section, and this section uses the software Matlab to numerically simulate equation (3) to verify the above conclusions.
e initial values of the system are (0.7,0.2,0.1), and the parameters are assigned as c � 2, η � 0.00714, β � 0.1, υ � 1, b 1 � 0.8, σ � 1/7, α 1 � 0.04, and Λ � 0.00334, where q and δ are two variables. In this paper, we will take some random values to simulate the change of S, E, and I ratio when R 0 < 1 and R 0 > 1, respectively. When q � 0.9, δ takes a random value between 0 and 1 so that R 0 < 1 holds. When δ � 1, q takes a random value between 0 and 1 so that R 0 < 1 holds. Its value of R 0 is less than 1 (R 0 < 1) when q and delta take the values in Table 1. Its value of R 0 is greater than 1 (R 0 < 1) when q and delta take the values in Table 2.
e global asymptotic stability of both disease-free and endemic equilibrium points can be seen from eorems 1-5. Figures 3-6 illustrates the correctness of the results obtained from the above theorem.
If R 0 < 1, by numerically simulating the stability of the disease-free equilibrium point, we can analytically conclude the following points. Figure 3 shows the variation curves of S, E, and I as the parameter δ varies when q is constant. e variation of S increases, and the rate of decrease and the rate of increase both increase with the increase of δ. e rate of decrease of E and I decreases with the increase of δ. Finally, S, E, and I all stabilized (see Figures 3(a)-3(c)). e change of S decreases, and the rate of decrease of E decreases with the increase of q. But the change of I has no obvious pattern with the increase of q. Finally, S, E, and I all tend to be stable (see Figures 4(a)-4(c)). Figures 3 and 4 illustrate that when R 0 < 1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R 0 < 1, the disease-free equilibrium point is tending to be stable. In Figures 5 and 6, the curve of S, E, and I variation with q and δ is simulated. By comparing the trend of S, E, and I with q in Figure 5 and the trend of S, E, and I with δ in Figure 6, we can draw the following conclusions. Increases in both q and δ lead to larger fluctuations in S and delayed stabilization times. It will lead to an increase in the value of S as it stabilizes. What is more, as δ increases, the rate of increase and decrease of E decreases. e monotonically increasing interval of I decreases and the rate of decrease of the monotonically decreasing interval decreases with the increase of δ. e maximum peak value of both E and I decreases with the increase of δ. And the time of reaching the maximum peak value of E is delayed and the time of reaching the maximum peak value of I is advanced. Furthermore, as q increases, the rate of increase and decrease of E decreases and E reaches its maximum peak earlier and at a lower value. As q changes, there is no obvious pattern of change in I, but it can be concluded that q has a significant effect on the change in the value of I. Figures 5 and 6 illustrate that when R 0 > 1, the changes of q and δ have significant effects on the trends of S, E, and I. And when R 0 > 1, the endemic equilibrium point is tending to be stable.

Model Simulation.
e real data of China in this paper are from the official website of the WHO. We preprocess the data, remove some bad points, and make the time continuous so that the data can be better adapted to the SEIR-QV model. e Chinese COVID-19 outbreak is fitted based on the SEIR-QV model. e parameters used in the model are listed as follows: c � 0.219, η � 0.00714, q � 0.5, δ � 0.2, β � 1 × 10 − 8 , υ � 0.275, b 1 � 0.05, σ � 1/7, α 1 � 0.035, Λ � 0.00334, b 2 � 0.05, α 2 � 0.035, and r � 0.9493. e simulation results and the actual data are shown in Figure 7.
e Chinese government has decisively taken strict preventive and control measures following the emergence of the COVID-19 outbreak.
ese prevention and control measures, along with vigorous advocacy by the Chinese government, have reduced exposure rates c and increased isolation rates q and minimized the development of aggregated outbreaks. e COVID-19 outbreak has become manageable and stable. As can be seen in Figure 7, the simulated data fit better with the real data from January 20, 2020 to May 20, 2021. However, there is a small outbreak of the epidemic in China in June 2021. e main reason is that the epidemic is out of control in many places abroad, and people at home are taking it lightly and weakening the efforts of prevention and control, thus giving the opportunity for the epidemic to spread.
As can be seen in Figure 8, both the isolation ratio and vaccine coverage play an important role in the control of the outbreak. e highest number of confirmed cases per day in China decreases as both q and δ increase. e correlation coefficient is used to measure the accuracy of the fit, and the value of R is found to be 0.9394. erefore, the results show that the SEIR-QV model can be used for COVID-19 epidemic development status assessment and has important implications.

Conclusion
An SEIR-QV model based on vaccination, isolation strategies, and the impact of different parameters on the development of infectious diseases is developed. e equilibrium point and stability of the new model are proved by using the basic regeneration number and the Lyapunov function theory. Simulation experiments show that the new method has certain theoretical value for analyzing and predicting the development of the COVID-19 epidemic. We can get the conclusion that there is a significant impact on the development of infectious diseases by different vaccination and isolation strategies.

Data Availability
e dataset used to support the findings of this study is available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.